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ABSTRACT 

This paper describes a microcomputer which was developed 
to compute ground radio beacon position locations using satel- 
lite measurements of Doppler frequency shift. Both tlie compu- 
tational algorithms and the mici'ocomputer hardware incorpo^ 
rating these algorithms are discussed. Results are presented 
where this microcomputer in conjunction with the NTMBUS-6 
Random Access Measurement System (RAMS) provides real-time 
calculation of beacon latitude and longitude. 
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SATELLITE DOPPLER DATA PROCESSING 
USING A MICROCOMPUTER 


1,0 INTRODUCTION 

A well established application of neax’ Earth orbiting satellites is the task 
of data collection and distribution, If the data is relayed in real-time the col- 
lection area is limited to the satellite field of view. For a nominal polar 
orbiting 1000 km altitude satellite such as NIMBUS-6, tliis corresponds to a 
maximum coverage circle of appi'oximately 6000 km in diameter. Since the or- 
bital period is on the order of 110 minutes worldwide coverage is obtained at least 
every 13 orbits. If the data is stored onboaixi by means of a tape recorder and 
"dumped” on command to the interrogating station, coverage is limited only by 
ground beacon deployment. NIMBUS-6 can operate in both onboard storage and 
real-time data relay modes, hi either case the data format seen by the ground 
processing computer is the same. 

Ever since launch (June 12, 1975) NIMBUS-6 has been relaying data from a 
number of platforms associated witli meteorological monitoring balloons and 
buoys (Ref. 1, 2, 3). In these applications platform location as well as platform 
data (for example, temperature and pressure) are required. The microcom- 
puter discussed in this paper has been evaluated by using both taped data as 
well as tying in directly to equipment capable of receiving and demodulating 
real-time NlMBUS-6 telemetry and RAMS Doppler frequency data. The teleme- 
try decommutation is rather straightforward and will not be discussed at length. 
The Doppler processing for purposes of beacon position determination is con- 
siderably more complex and it is this aspect of the microcomputer which will 
be discussed in detail. The same computer technology is directly applicable to 
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i'ortlicoming missions such as tlie near Polar orbiting TIROS-N Satellite Aided 
Seai*ch and Rescue demonstration. This TIROS-N demonstration, which is 
scheduled to begin in early 1981, will relay the Doppler shifted signal from low 
power "distress beacons" at 121.5, 243 and 406 mHz. In a manner tinalogous 
to tlie NIMBUS-6 buoy tracking, the ground based computer will combine the 
Doppler infornaation with known satellite orbit information to compute the 
latitude and longitude of the "distress beacon" (Ref. 4, 5, 6). 


2.0 SYSTEM DESCRIPTION 

The radio beacon-to-satelUte geometry is as indicated in Figure 1. As will 
be shown, the use of geocentric Earth fi.\ed coordinates greatly simplifies the 

computational problem. Bi such a coordinate system the range R, and its time 

derivative R can be expressed asj 




(X - Sj)2 > (Y - 3^2 (Z - S 3 V 


( 1 ) 


. . .^R . ■ Si) . (Y - S 3 ) ^ (Z - S 3 ) , 


In practice a constant bias terra must also be included in equation 2 to 
account for freq’iencj'^ offset. This offset is solved for along with tlie beacon 
position. The beacon coordinates Sj, Sv,, S^, are easily transformed to longitude 

and geodetic latitude. The satellite coordinates X, Y, Z and corresponding 

» » • 

velocity components X, Y, Z, being "Earth Fixed" inherently include tlie effect 
of Efirth rotation which is on the order of 7.3 x lO-^ radians/sec. 

In the NIlViBUS-6 Rams system tlie beacon is self timed and transmits at a 
nominal 401 . 2 MHz once per minute for a one second period. The one .second 
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transmission includes 320 milliseconds of unmodulated carrier followed by 
640 milliseconds of PSIC modulated digital data at a 100-bps data rate. For 
position location purposes NIRlBOS-6 measures the uplinlv Doppler shifted fre- 
quency associated with the unmodulated carrier* With reference again to 
Figure 1, the range rate, R, is related to the one way Doppler shifted signal by: 


where for NlMBUS-6 RAMS; 

fx - beacon tnmsmit frequency = 401.2 MHe 
R = beaeon-to- satellite range rate -7 lira/sec R 7 kni/sec 
C = speed of light - 3 x 10^ km/ sec 

As a conspquence of this Doppler shift and the frequency tolerance (±5 KHz) 
of the platform crystal oscillators, the NlMBlIS-6 receives a signal somewhere 
in a band of ±15 KHz, This frequency dispersion allows tlie RAMS system to 
track up to 200 beacons during any given satellite pass. Additionally each beacon 
is identified by a unique ID. The ID uses 12 bits of the 64 bits of data trans- 
mitted once per minute. Since the transmit frequency of the remote beacon is 
seldom known to better than a few KHz, the position location algoritlim must 
accurately determine tins frequency (or equivalently range rate) bias in order 
to compute beacon latitude and longitude. E.xpericnce has shown that tlie com- 
puter can recover this frequency bias to a high degree of accuracy (i.e. to better 
than a few Hz) even tliough data from the Doppler curve's point of inflection is 
e,\cluded. Figure 2 indicates a typical NIMBUS-6 Doppler pass. The recovered 
frequency bias in this particular case was on the order of 2 IxHz. The ordinate 
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TYPICAL ONE WAY DOPPLER MEASUREMENT 


12r* 



NOTE: FAIRBANKS STATION TO NIMBUS-6 
fo= 401.2 MHz 

START TIME 10 AUGUST 1975 
21 H 12 M 31,44 S 


Figure 2. Typical One Way Doppler Measurement 
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labeled "true zero Doppler" corresponds to the time when the beacon-to-satel- 
lite distance was a minimum. It should be reiterated, however, that as a con- 
sequence of the least squares algorithm discussed in this paper, Doppler data 
in the region of closest approach is not a requirement for frequency bias 
recovery. 

The overall system configuration used during the microcomputer evaluation 
is shown in Figure 3. The ground texmmal is a NASA prototype of a folly auto- 
mated NIMBUS-6/TIROS-N data collection station. For TIROS- N trackii^ the 
azimuth/elevation antenna drive will be directed by the same microcomputer 
which provides platform position location and telemetry data readout. 

3.0 COMPUTATIONAL ALGORITHMS 

The algorithms presented herein apply strictly to stationary platforms. 

Any platform motion will of course modify the Doppler signature and if this 
motion is left unmodeled a position location error will result. For slowly 
moving platforms such as drifting buoys this error is small and position deter- 
minations from successive passes can be used to estimate drift velocity. 

The error due to immodeled platform motion can be shown to be 
on the order of 0.2 km per km/hctur of uncorrected speed. At present the 
microcomputer is programmed to provide beacon latitude and longitude based 
on either a single pass or two successive satellite passes of Doppler data. For 
the single pass solution a minimum of 3 Doppler measurements are required 
since 3 unknowns (latitude, longitude and frequency bias) are solved for. Since 
the RAMS Doppler data rate is one measurement per minute, a maximum of 
approximately 15 data points per pass can be expected. Any additional informa- 
tion such as in situ observations of platform heading or speed, if sent back to 




Figure 3. 
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the microcomputer, can be used to construct a complete time history of position 
location over the pass interval, 

3.1 A Priori 

An "A Priori" or "first guess" of position location must be used as a 
starting point for any least squares type of algorithm. The better the initial 
estimate the fewer the number of iterations required to obtain a converged 
solution, in many cases such as that of a drifting buoy a previous position de- 
termination is the obvious choice for a first guess. However, there are cases 
when no prior beacon positions are available. One such application is satellite 
aided search and rescue where the beacon is associated with a downed aircraft. 

In this case, one straightforward approach which works well Is simply to 
let the computer scan a grid of the Earth surface in view of the satellite during 
the Doppler data collection interval. Each point on the grid is used successively 
as a starting value for the least squares algorithm until one that gives conver- 
gence is obtained, Another useful starting position Is talien to lie on a great 
circle on the Earth’s surface perpendicular to the ground track of the satellite 
at the midpoint of the pass. The point on the great circle Is selected at some 
arbitrary location between the ground track and the radio horizon. If data is 
collected during the time of closest approach (TCA) the following equations can 
be solved to obtain an estimate of beacon position: 

Sj^ Sj + (4) 

XSi + YSj ZS 3 = XX + YY ^ ZZ (5) 
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I 

•7 


( 6 ) 


XS, + 


YS2 + 


zs. 


- R2. 1 

min I 


A spherical Earth is used in this approximation with the geometry indicated 
in Figure 4. The satellite coordinates are Earth fixed geocentric and hence in- 
clude Earth rotation. The beacon-to- sate Hite range at time of closest approach 
can be approximated {Ref. 7) neglecting the effect of Earth rotation by; 


R 


min 



+ 


r 


R2 

V4 


+ P 



(7) 


where; 

a = Earth radius 
V = satellite speed 
h = satellite height 
P = a + h 


and 



f^i is the slope of Doppler curve at TCA {i.e. maximum slope of Doppler 
curve) 

Inspection of equations 4, 5, and 6 indicate the intersection of two planes {a 
straight line) in turn intersecting a sphere. The two single pass solutions thus 
obtained are approximations to the true and to the so called "image” solution. 

Ih tlie absence of Earth rotation the two solutions would be indistinguishable. 
Separation of these two solutions on the basis of Earth rotation effects is dis- 
cussed in the next section entitled "ambiguity resolution". 
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SATELLITE (X,Y,Z,^t,y,Z) 



Figiu’e 4. SimpUfied Solution Using Point of 
Closest Approaeli 
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Finally there iire other more accurate but more complex approximations 
which directly solve the range-rate observational equation independent of the 
relative beacon-to-satelUte geometry. Such schemes include the fitting of a 
high order polynomial to tl»e data or fitting portions of the data with segments 
of parabolas. However, experience to date indicates that oven a crude estimate, 
in error by several thousand kin, will for most cases result in a converged 
least squares solution in 5 iterations or less. The time for computing a typical 
10 point Doppler solution using 5 iterations on the present microcomputer is on 
the order of 3 minutes. 

3.2 Ambigtiity Resolution 

For any given satellite pass imd beacon-to-satelUte geometry there will be 
two locations, each at the same ortiioginal distance from tlie ground trace, where 
the riuige-rate (and hence Doppler shift) will be zero at TCA, While these 2 
Doppler chai'ae ter i Stic curves go tlmough zero at the same time, because of 
Earth rotation the shape of the two characteristics wdll be slightly different. 
Hence the measured Doppler data will be best fitted by the calculated Doppler 
corresponding to the side of the true location. On a sii^le satellite puss the true 
location ciui then be sepiu’atod from tile image location by comparing tlie RMS 
observed — minus-calculated residuals for each side. Simulations have shown 
that in tlie absence of data noise and in the presence of a well distributed data 
set this distinction is easily made. If two or more satellite passes of data are 
available this so called ambiguity is resolved simply from geometric considera- 
tions, Some insight into how Eiu'th rotation affects tlie Doppler characteristics 
can be obtained by writing equations in terms of beacon latitude, >/', and longi- 
tude, As a further simplification assume a spherical Earth or radius 'a* 


and a Polar orbiting satellite at height, h. The range-rate can then be ex- 
pressed (Ref. 8) as: 


R 


(a + h) a 
2R 


^ [c^s ■** t + ^] 


( 8 ) 


(OJ^ - <ij^) cos<j& sin [(o!g " <^0) t ” ^] “ 2c<jg sin0 


cos O)^ t 


where: 

= Earth rotation rate - 7.3 X 10 radians/sec 
oj‘ = Satellite angular rate = 10 radians/sec 
R = beacon-to- satellite slant range 

The effect of Earth rotation is nearly qrthoginal to satellite velocity at the 
equator and nearly coUinear with satellite velocity at the poles. The difference 
between the im^e and tnie Doppler characteristics is thus a minimum at the 
equator and maximized at the poles. Equation 8 ean he used to obtain an approxi 
mate analytic expression for this Doppler (range rate) eharaeteristio difference , 
6R, at the equatorwhere ,p = 0 and a circular polar orbit Is again assumed with 
X referenced to a TGA overhead sateLlite equator crossing at 0* longitude: 


2a (a + h) sin X 

5R - V - • - 

[a^ + (a+h)^ “ 2a(a + h) cos X cos t]^ 

where; 

;X.- ■■■: :'X =:beacon■l0ngitude^^■; 

t = time relative to TCA 
5H = true minus image characteristic 

and all other symbols are as previously defined 



As per equation (3) the frequency characteristic difference will depend on 
the beacon transmission frequency. A somewhat more general ^proach was 
taken via computer simulation and many computer plots of this effect for various 
satellite inclinations (i.e. including other than poiar)^ lati1udes of TCA, and 
beacon'-to-satellite TCA elevation angles have been obtained. Figures 5i and 
7 Indicate typical computer plots of the range rate difference due to Earth rota- 
tion for a satellite at 1000 km altitude» circular orbit and 80° incUnatloni andele-- 
vation angle at TCA of 60® at latitudes of 0® , 30® and 60®, It should be noted that 
at the equator the maximum excursion is on the order of 25 meters/sec (33 Hz 
at 401.2 MHz) and at latitude 60° it is 140 meters/sec (190 Hz at 401.2 MHz). 

The times indicated in Figures 5, 6, and 7 are relative to TCA and run to the 
radio horizon for the given satellite pass. Since the observed Doppler noise on 
the 401.2 MHz RAMS system is typically 0.5 Hz, Earth rotation signatures are 
observable. 

However , it should be pointed out that the degree of observability Is also a 
function of the quantity of data collected. For example 3 Doppler data points at 
a one per minute rate near the beginning of the pass may be insufficient to allow 
image discrimination. 

The real test, in the absence of a complete data span, is to generate a 
simulated "noisy" data set for one side of the pass and compute both solutions 
using the lea^ squares algoriihni and compare observed-minus calculated 
residuals. A number of such simulations have been rim and indications are that 
even under somewhat marginal condition of geometry and data quantity the true 
location can generally be separate from the ambiguous location on a single 
,pass.!-. ^ 
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SATtLLlTE ORBIT INCLINATION - 80 DEGREES 
UTITUDE OF BEACON -3D DEGREES 



-8“S'-4-20 2 A 6 0 

TIfC IN MINUTES 

Figure 6. Earth Rotation Effect 
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SATELLITE ORBIT INCLINATION - 80 DEOCES 

Latitude of beactn -go degrees 



- 8 - 6-4 "2 0 2 4 6 0 

TltC IN MINUTES 

Figure 7. Earth Rotation Effect 
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3.3 Atmospheric Effects 


As the signal from tlic beacon to the satellite traverses the Eiu’tli's atmo- 
sphere the Doppler frequency nK'asurement will be altered. The atmosphere, 
for purposes of radiowave propagation, cim be seijarated into two parts: the 
troposphere (sea level to 30 Um) mid the ionosphere (85 km to 1000 km). Tropo- 
sphere refraction effects are generally frequency independent up to about 30 Cdlz. 
Ionospheric refraction is frequency dependent varying inversely as the square 
of the operating frequency. The atmospheric bias on lX)ppler results from a 
scanning by the ray path tlirough media where the inde.x of refraction is other 
Uian unity. This atmo.spheric bias is a function of index of refraction along the 
path, elevation migie mid elevation migle rate, 

TROPOSPHERE 

If left uncorrected the maximum tropospheric induced rmige rate error is 
on tlic order of iO cm ^soc (0,13 Hz at -101.2 MHz) for a 1000 km altitude cireu- 
liu’ orbit satellite. This maximum effect will occur at low elevation angles for 
near overhead passes where elevation angle rates are highest. For systcjivs 
such as the HAMS where HMS frequency noise is on the order of 0.5 Hz, this 
biasing effect is small in comparison. However even a rather simple algorithm 


can remove most of this tropospheric bias at elevation angles tibove 10“ (Ref. D). 
An estimate of rmige rate bias cmi bo obtained using a series of well estjiblishod 
approximations (Refs, 10, 11). In the troposphere the group mid phase velocities 
iu*e equal in magnitude mid of the same sign. The range rate bias estimate can 
thus be obtained by' differentiating the expression for rmige bias with respect to 
time. If an exponential atmosphere is assimied the residtmit rmige rate bias 
estimates cmi be quickly derived. That is — 
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( 10 ) 


AR ■ -—j N, 


e"»*’'dh 


10-® 
k sin E 


Where: 

- 1 + Nj. (10“*') surface index of refraction 
Ns - 350 

K - line-of-sight elevation angle beacon-to-satellite 
k - decay factor 1.0 x 10-* meters"* 

:uiid 


NR. 


'WR ^ 

Nt 


^'os 

ksiu'E 



( 11 ) 


Now the elevation rate corresponding to an overhead circidai* pass Ci\n be 
written as: 

±tjj. 

E — (12) 

(, . ' ..iiit 

\ 11= sin= E ) 

where ^ ' 

E elevation rate (positive sign for Increasing elevation) 

h satellite height - 1000 km 
a - CiU'tli radius - G378 km 
-'s satellite tuiguliu* rate 10"^ radums/sec 
For the values indicated the elevation angle rate is 1,4 (10 ~^) radians/sec at 
10*’ elevation iind as previously stated using equation 11 the range rate bias is 
seen to be approximately -10 em/sec. As the elevation angle approaches 90“ 
(i.e. directly overhead) tlie troposphere range rate bias tends to zero. 
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IONOSPHERE 


A siinilai' approxinuUioii cjin also be derived for liie ionospheric liopplcr 
effect (Ref. 12). For frequencies tibove about 100 MHz the inde.x of refraction 
within tlie ionosphere can bo expressed as; 


N.. e - \ ■ 


I - 


‘ .1 m 


0 


40. 3 N,,' 
i'i 


(13) 


where: 

Electron density (electrons ■‘'meter'*) 

e Electron charge - l.li02 x 10"*‘’ coulombs 

m Electron mass i).ll \ 10''” kilognuns 

t j, Free space dielectric constai^t 
* 1 ' . 

8.S55 X 10 ' tarads ' meter 

2 ' f (radimis sec) 

It should bo noted that witliin the ionosphere the index of refraction is less 
tlum unity tmd during a given satellite pass tliis results in a Doppler bias of 


opposite sign relative to tropospheric bias, 'I'lie following expression includes 
a correction for horizontal gradient effects (e,g. ilay-to-night trmisilions) if 


information regarding the ehsmge in vertical oleelron content. 1^., is avaihible in 
terms of the local elevation angle, E*, at tlie height of electron density nuixi- 


mum, h,„. 


HOKIZONTAL 

fiRADlEET 


M0.3 1 . cos 


VR. 


t ■ sin 


1 - 



meters see 


(M) 


It) 



where: 


f LOOO KM 

= I Np dh (Electrons/meter^) 

8 5 KM 

f - frequency of transmission (Hz) 

E* = local elevation angle of Une-of-sight at height of maximum density 
E* = local elevation angle rate of line-of*- sight at height of maximum 
density (radians/see) 
and from geometric considerations 


E* 


arccos 


a + h 


cos E 


(15) 


where: 

a = Earth radius (or spheroid radius of curvature if available — 
appendix A2.2) 

h^ - height of maximum electron density = 350 km 
E = elevation angle at beacon 
E* = local elevation angle at h^ 


E+ = 



(16) 


where: 

E = elevation angle rate as observed at beacon 
E* - elevation angle rate at h^, 
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The vertical integrated electron content can vary between 10^® and 10^® 
electrons/meter ^ depending on geographical location, season, time of day, sun 
spot cycle and so on. The continuously updated NASA world ionospheric model 
developed for providing ionospheric corrections to NASA tracking data is one 
source of integrated content values for offline corrections to beacon data. The 
expression given by equation 14 can be used to provide estimates of maximum 
ionospheric effects. For purposes of comparison with tropospheric effects a 
ground elevation, E, of 10° is again assumed along with a ground elevation rate, 
E, of 1. 40 X 10'^ radians/sec, A typical daytime vertical integrated content of 3 x 10 * ^ 
electrons/ meter ^ with no horizontal gradient is also assumed. Using the fore- 
going, this results in a maximum 401,2 MHz range-rate bias of 35 cm/sec 
(Doppler bias -0. 5 Hz), At 100 MHz this bias increases to 6 meters/sec 
(floppier bias -2 Hz). 

Since system frequency measurement resolution will typically be 1 Hz or 
better, an ionospheric correction in the 100 MHz region is imperative. At 
400 MHz and above, such a correction is optional but should be included if full 
system capability is to be realized. While the atmospheric effect has been re- 
ferred to as a bias, it actually has a time varying characteristic over the full 
satellite pass. It therefore cannot be as easily solved for as, for example, the 
oscillator frequency offset which is nearly constant over the satellite pass. 

Using beacon frequencies higher than 400 MHz will further reduce iono- 
spheric Doppler effects. However, if simple low gain (and hence broadbeam) 
antennas are used at both satellite and beacon to optimize coverage (Bef. 13), a 
nominal link frequency of 400 MHz is probably a good choice. 
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3.4 Orbit Update Procedure 

Accurate recovery of the unknown beacon's coordinates is directly dependent 
on the accuracy of the reference orbit. In practice it is found that the absolute 
accuracy of the reference orbit deteriorates with time resulting in a corre- 
sponding degradation in the accuracy of the recovered beacon coordinates. The 
most significant error In the reference orbit is the one associated with the 
’’along track" error. This error may be thought of as an error in the mean 
anomaly or alternatively, as an error in the epoch time of the reference orbit. 
Using Doppler daty from a beacon of known position the orbit update procedure 
solves for this "alo.ig track" error and corrects the reference orbits epoch 
time accordingly. Each range rate (Doppler) observation is considered to be a 
function of the reference orbit's epoch time T, and the range rate bias B. The 
linearized observational equation is then given by, 

where , 

R° = observed range rate 

R = "computed" range rate — based on reference epoch T, and beacon 
position 

^ = partial derivative of "computed" range rate with respect to epoch 
time T- 

\T ~ correction to epoch time T 

AB = correction to range rate bias 

As a practical matter the partial derivative is computed numerically as a 
first difference, — that is. 
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( 18 ) 


R(T + AT^ - R(T) 

' AT 


Experience has shown that for the geometries considered in this paper, a value 
of AT equal to one second is a satisfactory interval when using the approxima- 
tions indicated by equations 17 and 18. The quantities AT and AB are then 
solved for using the method of least squares. Each data point is processed for 
inclusion in a set of normal equations designed to yield the updated orbit solu- 
tion. Let N represent the 2 by 2 normal matrix with elements N. . , The four 
matrix elements are evaluated (e.g. Ref. 24) as follows: 





similarly the 2 by 1 matrix 



is evaluated as: 
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Again all suma aiv carrioil over all n data points yintv a scparato linoari/.otl 
equation exists at eneh measvm>mont time. 

I'he least squares aotutiun is obtained by solving this nortnal system of 
linear equations mui is given by, 


rAT'i 


\u 


N-* C 


luul linally. T i;, * \T juid H It, * \H 

represent the e\irrent solution values (Tj, and be»ng initial values) wlueh are 
useil as starting values for the next iteration, 'I'his iterative process is eon- 
si<lered to liave eonverge^.i when the absolute value of the \'V and \B eorreetions 
are aeeopttibly small. 

I'his newly solved for value td’ orbit epoch lime T, is then used to update the 
reference orbit on all subsequent solutions for unknown beacon local ions, 
ilefinite impri'vetueul in the recovery of beai’on coordinates results if tl\is prt>- 
t'eduiv is used on each t>rbil:il pass. This requires that al least one reference 
beacon be avaitat)le mi each puss. The procedure itself is automatic, the only 
information needed is the identification number of the surveyetl in refereiu’c 
beacon which as a matter of couveuioiu'e mm be loeated at tlie site receiving the 
satellite relayed Boppler data. 
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rho noouriu'y ol‘ boacon pi^stliou location compntaliona will bo innuonccil by 
IX>pplor data iiualily ami tpuuitlty; boacon-to- satellite geometry; ami accuracy of 
satellite ephemeris. Mcaes for upgrading the salellite ephetueria ^Yere dis- 
cussed in the pri'v ious section, Kxcept for the near overhead pass case, aatel“ 
liti' pt>sition errors eou|ile approximately into equivalent beactm yiosition errors. 
That is if file orbit is well ileterminei.1 (known within tens of meters) the contri- 
bution tci i>osition location error will be in the tens of meters. Krror studies 
(c.g. Kef. M) indicate, howevi'r. for the near overhead case ptisiiion location 
errors become as large as ten limes the satellite position error. 

lloppler liata processing experience has shown that for the single pass case, 
Lu'ju'on height above.' the gcoiil is generally not observable. It ejm be shown that 
this error translatos into a position I'rror approximaloly equal to the height 
unce'rlaittly multiplied by tan K, where K is the beaei>u-to- satellite elevation 
angle. If more’ than one' pass of elala is avail.'U.)le fntm a slatieinary beacon, 
hoiglu e’an bo sedved for along with latitude. Kmgitiule and frequency bias. ’I'he 
nuuinunn uumbe'r ed' data p(>ints re'quire'd bn* a pe'sitiem detormiuatie'n will of 
cemrse eqvial the number ed' uuktu'wns being solvi'el fen'. 


d1ie data quality is generneil by siu'h faele)i-s as time lag aeeuraey and 
fre'queney stability Lhianiglunit Uk' system. 'I'ime tag errors eeniple in d'rectly 
and are inelistingiiishable from along track e>rbil ermrs. For example, a 


0,1 second lime error is equivalent in elTeet U) mi 700 meter along tr.ack orror 




since the spacecnil’t velocity is appi'oximately 7 km see. 'I'he MlMlil'iS-ll 
TWEKLK^llAMS system was designcLl for a if> km sint;lo pass position accuracy 
imd hence a time tay: accuracy of lO.l seconds was adequate. The TUK>S-N 
lime tag for im equivalent Doppler tracking system will be lO.OlO seconds. 

Another importemt factor is the i>verall frequency sl:ibnity including beacon 
trmismitter and all hetrodyning oscillators up through the average frequency 
measurement, Clenerally the HAMS system provides sufficient data (c,g. 10 or 
nu>re Doppler points at a one per minute rale) such that the noise (typically 
0.5 Uz at '101.2 MUz) mul the solved for bias i>ffscl (typically several Kllz) arc 
of little or no consequence in terms of location accuracy. The noise is avcragcil 
by the least squares fitting technique and the bias is easily e.xliacted U> wilhin 
system resolution values (i 1 llz). However long term oscillator ilrift. which 
chtmgcs the slope of the IX^ppler curve, will introduce position errors. Doth 
simulations and error tmalyses indicfite that the HAMS sensitivity to frequency 
drift for a full Doppler data span can be estimated by a simple e.xpression de- 
duced from an approximate differenliation of equation (7) aamoly: 
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where; 
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- ixjsitlon error (meters) 


eos E 

Alt beacon- to-salellite error (meters) 
c speed of light (meters ''see) 
f bvaeon nominal frequeney (Hz) 

Af frequency tlrift (llz sec) 


E 


b eae o n- to- s ate 1 1 i to e I e v at i o n 



V .satoUiio spivii 7.1^5 'lo' motors soo 
k iyploiil miiMinnuu iU’ooloration 25 motors ' sec 
Uv obsovving ohaut^os in troquonoy bins on eonsooutlvo pnssos, Af Inis boon 


obsorvotl lo bo tyv>iouUy i li/, lU'r 1000 soooiuis lor tho HAMS system. The ox- 
pooioO |>osiUon error oorrospomling to such a ilrift at 401,2 MU/, for a TCA olo- 


vjitlon of Ifi'’ nill tivoroforo bo on the order of 100 motors. 


At a TCA elevation 


ot so" this error wonUl Ineroaso to 400 motors. 

rho otToet of tiu' aUnosphoro is to ohango tho slope of tho Hopple r ourvo and 
in this sons<' introtluees an error in |x>sitlon loeatton very similar to that due to 
frequeney drift, I'or example, as indieatod in the tnwlous soetion, an iono- 
spherie exeursion of u2 Ih* at 100 MUv* oan easily bo expootetl for a near over- 

ho;ul pass. If this ehange oeeurs over a typieal 1000 seooml puss an "equiva- 
lent ilrift" of 0.004 11/ see is exv'torionood. Putting this value into oquaUou 21 
and as.^undng a Tt'.A elevation of StV’ indtoalos an oxpeotod posit ion error of 0 km. 

At 401.2 Ml!/ this effect for an uneorreoted daytime ionosphere ("oqulxqilenf 
drift" of 0.00! 'seel and the same geometry would be on tho order of 0.4 km. 
If even a nondnal oorreetlou is made based upon equation 14 this error can be 

reduc'd to 20'\‘ of the imeori'eefed value, two frequency selu'me such as 
associated with the Navy Transit system can be employed to provide real-time 
ionospheric eorreetions which even during high solar activity reduce this error 
to a value on tho order of 10 meters (Uef. 15). If only a portion of the Doppler 
curve is used for position determination the error due to frequeney drift is re- 
duced. However in this case the position em*or due to data noise increases to 

the extent find frequeney noise often becomes the predominant system error 


source. 



-1.0 MK'lU>(.'OMPl! I'KK I^KSCHIPTIOK 


Evoi* siiifL' tlio 107r> NlMlUlS-0 Imuu’h, tho t’ompulatioiis ri'qiiirL'd for posi- 
tion location of Lnun’s, mctoorolo};ic:il balloous, and other experimental pla1h>rms 
uain^’ the NIMlirS-tl system have been aueeessfully carried out using a 

centrally located eonvenlioniil general purpt'se computer. Doppler position 
location systems similar to the NlMlU'S-0 scheme will be carried by one or 


more of the forthcoming series of TIllOS-N spacecraft. 


The TIKOS-N platform 


monitoring :uid pi)Sition loeatitm system will be applied to btHit user eommunil> 
platforms and "distress beacon" Uaailions during the Satellite .\iileii Search and 
Ueseue Deimmstrution. Much interest has been expressed by the NASA usei' 
community' iit being able to set up low cost terminals capsible of receiving and 
processing NllillU’S-i! and TlllOS-N telemetry and noi>pler data direetl\ at the 
user's >ite. U is in direct response to this user need that the prototype miero- 
eomputer doseribed in this paper was developed. The required Unv eosl, eom- 
paet, anil reliable eomputer facility suitable for a self contained user terminal 


has been mailo [Possible through the advent of largo-seate integration (l.Sl) 
mie roeleelronie eii'ciiits tHef. U>) where thousands of circuit elements are in- 
tegrated in one ehii>, f'or example typical Hatidom-Access Memory (HAIM) chips 
measuring ;> by f> milliinclcrs provide an equivalent of o\er lii.bOO tnmsistors 
whei'e each "transistor cell" is cainibte of storing one binary bit of information 
(Hof. 17). 

In such chips the lime required to write one bit in an arbitrary 


location or to read it out is about 20U lumoseeonds. 


The remainder of this sec- 


tion will discuss how such LSI chips were intereonneeted ;md programnuni to 
perf-irm the specific task of eomiuttiug beacon latitiule and longitude using 



satellite observed Doppler frequency data and the known satellite orbit as input. 
Microcomputer alphiuuimeric output has been obtained using either a standard 
teletype printer or a catliode ray tube readout. The microcomputer input has 
been successfully interfaced and accepted data in a variety of ways including: 
paper tape via the teletype terminal; teletype keyboard; American Standai*d Code 
for hiformatiou Interchange (ASCII) keyboard for direct manual entry; cas- 
sette tape recorder output; and direct real time tie in to a NIMBUS-6 telemetry 
receiver. The next few sections will discuss tlie Central Processing Unit, 
Memory, and Input-Output Peripherals of the N1MBUS-6/T1ROS-N position loca- 
tion microcomputer. The total power requirement for this microcomputer is 
on tlie order of 7 watts (+5 V at 1.2 A mad +12 V at 0,1 A). 

The basic block diagram of the microcomputer imd interface is shown in 
Figure S. 

4.1 Central Processing Unit 

Every computer requires some type of central processing unit (CPU) where 
the actual computing is carried out. The CPU usually controls all tlie operations 
of the computer. Compaiters designed as recently as ton years ago used CPUs 
consisting of rows of printed circuit boards requiring a fairly large enclosure 
to house all CPU related components. Today tlie entire CPU can reside on a 
single LSI microprocessor chip (Ref. 18), At present tliere are over 50 different 
microprocessor chips on the market (Ref. 10). In tlie position location micro- 
computer the pai*ticulai* LSI chip incorporated is tlie 40 pin MOS Tocluiology, 

Inc. ty-pe MCS 6502 which incorporates an 8 bit data bus (Ref. 20). This chip 
provides a 16 pin address bus and tiierofore is capable of addressing up to 2‘® 
or 65,536 bytes of memory where a byte by definition consists of 8 bits. The 
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Figure 8. Basic Position Location Microcomputer and Interface 











term "bit’* refers to imything tiiat can be assigned a binary value of 0 or 1. The 
MCS0502 has the basic timing oscillator and clock driver on chip, thus elimi- 
nating tlm need for im external high level two-phase clock generator. This CPU 
chip cjm be used with jm externally generated time base consisting of either a 
transistor-ti’ansistor logic <TTL) level single-phase clock, crystal oscillator, 
or RC network. In tlUs application an external time base consisting of a 1 MHz 
crystal oscillator is used to generate tlie required two phase clock pulses. 
Nominal clock pulse wldtli is 430 nsec. As in most microcomputer systems, 
tlie timing of all data transfers is controlled by the system clock. The clock 
itself is actually two non- over lapping squiu'e waves. This "two-phase** clock 
system essentially provides two consecutive positive going switching pulses 
where the address Hues and road-write commtmds are set up during the first 
pulse mid data is trmisferred (cither read or write) during tlie second pulse. 

There are also two interrupt input lines to tiie microprocessor. The pri- 
mary purpose of tliese lines is to permit a mode of operation which minimizes 
the computer time spent in interrogating input porlpheral devices imd reading 
out via output peripheral devices except during times of specific interest. 

The internal organization of the processor can be split into two sections, 
hi general, the instructions obtained from program memory are executed by 


implementing a series of data transfers in the "register section’* of the chip. 


The control lines which actuaily cause tlie data transfers to take place are 


generated in the "control section". 

Instructions enter the processor on the data bus, are latched into tlie in- 
struction register, and tu*e then decoded along with timing signals to generate 


Uie register control signals. 
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The timing control unit keeps track of the specific cycle being executed. 

Each data transfer which takes place in the register section is caused by 
decoding the contents of both the instruction register and the timing counter. 
Additional control lines which affect the execution of the instructions are derived 
from the hiterrupt logic and from the Processor Status register. The Interrupt 
logic controls the ^>rocessor interface to the interrupt inputs to assure proper 
timing* enabling, sequencing* etc, which the processor recognizes and services. 

The Processor Status register contains a set of latches which serve to con- 
trol certain aspects of the processor operation, to indicate the results of pro- 
cessor arithmetic and logic operations* and to Indicate the status of data either 
generated by the processor or transferred into the processor from outside. 

The real work in the microprocessor is carried out in the register section 
of the chip. At 1 MHz* the data coming into the processor from the program 
memory, the data memory, or from peripheral devices, appears on the data bus 
during the last 100 nanoseconds of the second (or phase two) positive clock pulse. 
No attempt is made to actually operate on the data during this short period. In- 
stead, it is simply transferred into the input data latch for use during the next 
cycle. The data latch serves to trap the data on the data bus during each Phase 
Two pulse. It can then be transferred onto one of the internal busses and from 
there into one of the internal registers. For example, data beii^ transferred 
from memory into the ’’accumulator" will be placed on the internal data bus and 
will then be transferred from the internal data bus into the accumulator. If an 
arithmetic or logic operation is to be performed using the data from memory 
and the contents of the accumulator, data in the input data latch will be trans- 
ferred onto the internal data bus as before. From there it will be transferred 
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into the "arithmetic logic unit" (ALU). At the same time the contents of the 
accumulator will be transferred onto a bus in the register section and from 
there into the second input to the ALU. The results of the arithmetic or logic 
operation will be transferred back to the accumulator on the next cycle by trans- 
ferring first onto the bus and then into the accumulator. All of these data trans- 
fers take place during the Phase One clock pulse. There is also a "program 
counter" which provides the addresses which sequentially step the processor 
through instructions in the program which resides in memory. 

The accumulator is a general purpose 8-bit register which sorts the results 
of most arithmetic and logic operations. Ih addition, the accumulator usually 
contains one of the two data words used in these operations. 

All logic and arithmetic operations take place In the ALU. However, the 
ALU cannot store data for more than one cycle. If data is placed on the inputs 
to the ALU at the beginning of one cycle, the result is always gated into one of 
the storage registers or to external memory during the next cycle. Each bit of 
the ALU has two inputs. These inputs can be tied to various internal busses or 
to a logic zero; the ALU then generates the SUM, AND, OR, etc. function using 
the data on the two Inputs. The "stack pointer" and the two "index registers" 
each consist of 8 simple latches. These registers store data which is to be used 
in calculating addresses in data memory. The address bus buffers consist of a 
set of latches and TTL compatible drivers. These latches store the addresses 
which are used in accessing the peripheral devices (ROM, RAM, and INPUT/ 
OUTPUT). 

The microprocessor accesses position of memory for program instructions 
and f 1 to perform the overall computer Action, 
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4.2 Data and Program Memory 
DATA 

The memory section of the microcomputer is in two distinct segments. The 
data which consists of the Doppler measurements performed by the satellite as 
it passed over the radio beacon is stored in 4 k bytes of Bandom Access Memory 
(BAM). The term '^random access" means that any word in the memory may be 
accessed without having to go through all the other words to get to it. The only 
drawback of BAM memory is that when power is lost from the chip all memory 
content is lost. 

The Doppler data is loaded into the computer in any of a number of ways 
including: punched paper tape via the teletype, magnetic tape, real time NlMBUS-6 
receiver tie in or manual entry via the keyboard. Also stored in this memory 
bank are the coordinates and velocity of the satellite at some epoch time near 
the time of data collection. As indicated in section 3.4 the microcomputer has 
the capability of refining this orbit if Doppler data from one or more reference 
beacons is available. The satellite position and velocity over the Doppler pass 
are computed at the time of each data point and the memory need only store the 
6 numbers representing the satellite state vector at one point in time, hi addi- 
tion to this "data storage" the BAM also provides a certain amount of "workii^ 
storage". This is required for such things as storage of intermediate results in 
arithmetic operations and peripheral output data storage. The BAM of the posi- 
tion location microcomputer discussed in this paper consists of 32 t}rpe .^102 
RAM chips, where each chip has 1024 bits of memory. 
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PROGRAM 


Program algorithms are those discussed in the body of the paper and the 
appendix. Originally as the position location least squares program evolved, 
it was loaded into RAM using a cassette type magnetic tape. After the program 
design was achieved it was loaded into Read Only Memory (ROM) chips or more 
correctly into an Erasable Programmable Read Only Memory (EPROM). These 
particular chips are permanently programmed by means of a "PROM Programmer 
Unit". The memory can only be erased after a several minute exposure under a 
special high intetaity ultraviolet light. Once erased, the chips can be reprogram- 
med. This particular program memory consists of eight type 2708 EPROMs 
where each EPROM chip contains 1024 bytes of storage for a total program 
storage of 8192 bytes. 

It should be mentioned that a complete floating point math pack of subrou- 
tines is included in the EPROM to allow evaluations of all arithmetic algorithms 
including multiplication, division, trigonometric and exponential functions. 

The overall floating point microcomputer precision is 12 decimal digits which 
is more than adequate for this particular application. 

4.3 Input-Output Peripherals 

As already suggested a wide variety of peripherals can be interfaced with 
the microcomputer. Figure 9 shows the microcomputer interfaced with a key- 
board and CRT video display. Figure 10 shows a close up of the printed circuit 
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Figure 9. Position Location Microcomputer 



Figure 10. Microcomputer Printed Circuit Cards 
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cards ^^4llch make up the microcomputer along with the required d.c. power 
supply. One important consequence of this position location microcomputer 
development is a demonstration of the relative ease with which such a computer 
can be interfaced with the "outside world". The controlled transfer of data and 
program information between the processor and numerous peripheral devices 
(i.e. "handshaking") has been carried out. In fact, a Cross Assembler was 
used where output from a DEC PDP 11-40 directly programmed the EPROM in 
the microcomputer. Whenever changes were to be evaluated, a new EPROM 
was plugged in and programmed. When the program was considered obsolete 
the EPROM was erased using ultraviolet light and hence made ready for future 
programming. Of course, certain portions of the programming such as the 
"math pack", once checked out required no further modification. The computa- 
tion time for a typical Doppler pass of 15 data points over a 10 minute observa- 
tion interval is typically 20 seconds per least squares iteration. Usually 3 
iterations are sufficient to provide a converged or final solution. The computer 
recognizes the platform identification I-D (as contained in the telemetry) of each 
platform and sequentially stores the data in RAM as it comes. As many as 50 
platforms may be tracked at any given time by NIMBUS-6. Once the data is 
collected the operator "punches up" the I-D of those platforms he is interested 
in computing positions for. The microcomputer then scans the data portion of 
the RAM for one particular 1-D and at each data time computes the satellite 
position and velocity as well as the partial derivatives required for the least 
squares solution (see appendix). Once all the partials and satellite vectors have 
been computed for a given beacon the solution is obtained and printed out in 
alphanumerics either on teletype or CRT. Then the computer goes on to the 
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next I-D and so on. Experience has shown that given the same Doppler data set 
and oii)lt information the microcomputer eomputea the same (to 8 significant 
figures) position as much larger general purpose machines such as the IBM 360- 
95 or the DEC PDPll-40. Finally, if the input data is stored on punched paper 
tape, an optical reader has been shown to be an ideal low cost way to rapidly 
read in information to the microcomputer RAM. In this mode tlie data is read 
in as fast as the tape moves under an illuminating lamp. 

5.0 NIMBUS-6 RESULTS 

Numerous RAMS type beacon location computations have been performed 
using tile position location microcomputer. Also along track corrections to the 
nominal NIMBUS-6 orbit have been made using the reference 401.2 MHz beacon 
located at GSFC, Greenbelt, Maryland. The NIMBUS-6 spacecraft is tracked by 
the NASA Minitrack interferometer tracking network. The "definitive orbits" 
using this angle tracking data are provided about two weeks after the data is 
talcen and ai*e accurate to better tlian a few hundred meters. Long term (one 
montli or more ) orbit predictions for antenna pointing and data acquisition pur- 
poses will be significantly less accurate, hi either case real time orbit im- 
provements can be realized using one or more reference beacons, hi fact, if 
several such reference beacons are deployed, a real time geometric solution 
can be made and the need for orbit data is obviated. 

Figure 11 presents some typical beacon location results using the NIMBUS-6 
RAMS system and position location microcomputer. These particular passes 
were handled in real time — that is, as soon as the data was recorded the posi- 
tion location computations were performed. It has been observed that the 
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Figure 11. Typical Position Location Results 



iHipptoi' JjUsv nol«o vk'tcnuliiiHl tho UMS v'oiiUUuils of obeoi'vcti minus 
oomputud nppuudlx) Is typlotvUy loss Hum t Hr, at 401.2 MUa. This 

Is ounslslmvt with th« HAMS frequency i*csolvitlon of t l U?, corrcspmviUng iu a 
vangc-ratc resolution of il.5 motors,''soc. 


0,0 t'ONCUiSlONS 

A low cost, accurate and reliable special purpose n^lcrocomputcr for the 
real tln\e etnnpulatloxi of position locations of radio beacons ob.se wed by satel- 
lite has bei'u successfully develov>ed and testcxl. This basic u\icrocv>n\puter 
desijtu will be Incorporated In tNw NASA NlMllUS-tj '4'IHOS-N low cost Local 
User Terminals (LUT) currenUy belint aasetubled for demonstnvUon purposes. 
Vho LU^L will be totally automatic with the arlmuth^elevatlon aivtenua drive In- 
formation .also providwl by the ndcrocompvtter. The \)OSUlon locution of beacons 

A 

will avitomatlcally be pitudvlod oti a single p.ass and two consecutive pass basis, 
This allows for imprvned accuracy of fixed be.acous auvl v'osiUou civmge com- 
putallou of drifting beacvms (i.o. buv>ys). I'ho position eounnilatiotv accuracy Is 
dependent on orbit accuraey and l>i>ppler data quality tutdev various goometrle.s. 
Tho udcrocomputer computational accuracy la couslsient witl\ that t^bialned with 
much larger general purpose computers to vH vlecitual digits In floating point 
arithmetic. 
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Al’PKNDlX 

i.i;ast-squau microcx^mputer ai-gorithm 

The i'ol lowing relationships are programmed aloniv with the 'hnath package’' 
in sectors of Um microcomputer KPROM mid are calletl up as vequlroti. The 
particular mathomatlcal relation called up will depend on tlie type input being 
providixi. For example ii the satellite orbit int'cu'inatlon is already given in 
gt'oeentric Farth OxchI eoonlinates the E'arth rotation matrix computation will 
not be used and so on. The desired beacon position is always in Eai'Ui fixed 
cooi\linates (i.e. latitude mid longitude) and once the satellite position and 
velocity arc expressed in Earth fixed coordinates the time varying range rate 
(t)oviplcr) between beacon mul satellite is easily expressed (o.g. equation 2), 
Conveniitmal coordinate transformations (Ref. 21) are Incorporated in the 
niicrocomputer. A least squares (Refs. 22, 23, 24) technique is applied in tile 
data reduction to assure maximum position location accuracy. 
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A,1 COOKDlNA*i'K SYSTEMS 


A, 1.1 Geoeontric inertial of Date (CartesUm) 

Tilts system has its origin at Hie eentev of the Earth. The u luui v axes He 
in the equatorial plmie with tlie u iixls towards the vernal equinax of date. The 
w axis is positive timnigh the north pole, mid the v axis in a dirootion such that 
a right handed cooixHiuite system i i formed. 


A.l.ia Geoeentrie EarUi Fixed (Cartes Imi) 

The origin of this system is also at the center of tlie Eartli. As in Uie 
Geocentric hicrtial system Uie x and y axes lie in Uie equatorial pHme however 
the X axis passes Hirmigh tlie Greenwich meridian. The po.siitive a direction is 
again dirt^eted towards the north pule and the direction of Uie y axis is con- 
sistent vvlUi a right luuided coordinate system. 


A.2 COOKDINATE 'niAMSFORMATlOES 

.V.2.1 Geocentric Inertial to Geocentric Earth Fixed 
The transformation is given by* 


X 


” u 

y 

A(:'l 

V 

z 




Whei‘e is the angle between the line formed fixrni the center of the Earth to 
the true vernal equinox 4vnd Uie Geoeentrie Earth fixed x axis, 
is usually available as; 


* : (t to^ 
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and is the angle at 
to ~ 0 hours GMT of tlie epoch day and 


t - epoch time GMT 
e - Earth rotation rate 
The Earth fixed velocities are given by: 



in all of tlie above: 


cos sin c> 0 

A(i^) - “ sin cx>st* 0 

_ 0 0 1 


A-3 





A.2.2 Geodetic Latitude and Longitude to Geocentric Earth Fixed 
The Eartli shape and beacon location parameters are as follows: 
a - equatorial radius of Eai'tli reference spheroid 
e - reference spheroid eccentricity 

- geodetic latitude of beacon (positive north) 

\ - longitude of beacon (positive east) 
h ■ height of beacon above reference sphei'oid 
Let 1 represent tlie magnitude of the radius of curvature in tile plane of tlie 
prime vertical then: 

a 

(1 “ sill' v") ‘ 

:uid the beacon's geocentric coordinates iire given by: 

Sj (: t h) cos v' cos \ 


Sj - (;. » h) cos sin \ A“5 

S3 “ [: (1 - + bj sin 


A. 2, 3 Geocentric Earth fixed to Local Elevation 

The atmospheric corrections (section 3.3) require computation of local 
elevation, E, and local elevation angle rate, E. E and E are obtained as follows 
for tile local topocentric E:irtli fixed cai'tesian coordinate system 1 ’^; 

E - tan~* A-(> 


E 
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where x, y, z and Sj, S 2 , S 3 * are Geocentric Earth fixed coordinates corre- 
sponding to the satellite and beacon respectively. And the rotation matrix M is 
given by: 


M - 


- sin \ 

■ sin 4> cos 
cos 4> cos \ 


cos iV 

• sin 4> sin \ 
cos sin N. 


0 ~ 
cos 
sin 


also since the beacon is Earth fixed: 




r— . »- 
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which permits E to be computed as: 


E - 


(e2+ {-2)^3- 


'A 


+ qf 
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A- 11 


where the slant range, B, is computed from 

R2 = (X - Sj)^ + (y - + (z - 83)2 A-12 

The atmospheric corrections are not applied until after the first least squares 
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iterative beacon solution for Sj, 83. S3, (or equivalently o, V) to assure proper 
convergence. 
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A.3 LINEAR LEAST SQUARES SOLUTION FOR BEACON LONGITUDE 


AND GEODETIC LATITUDE 

The least squares approach incorporated in the microcomputer is conven- 
tional in tliat the equation of observation (equation 2 of text plus a constant bias 
term) is linearized. The result is a set of linear equations equal to tlte number 
of Doppler data points, m general the number of points (hence, equations) will 
exceed the number of unkno^vns which in this case is three (latitude, longitude 
and range rate bias). The procedure outlined in tl>e following indicates how tliis 
set of lineal’ equations are reduced to a new set, namely; the normal equations, 
where the number of equations equals the number of unknowns. The solution of 
the normal equations in the microcomputer is achieved in a predetermined 
number of iterations (usually 3). The resultant least squares solution is con- 
sidei’ed to have converged when successive iterations alter tlie result by im 
amount less tlian that attributable to data noise. 

By expanding tlie equation of observations in a Taylor's series about the 
a priori position and neglecting higher order terms; 




A 3R :^R 

^ ^ (v" - v-'o) ' ^ ~ ^o) + Cf* ■ ^o) A-13 


w’here 

R« ^ the observed range rate (Doppler) 

R - computed range rate baaed on current estimate (a priori) of beacon 
(see section 3.1) 

= initial geodetic latitude 
\q « initial longitude 
B(j = initial estimate of range rate bias 
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On the first iteration values of and are generally available to within a 
few hundred km based on prior knowledge obtained in such ways as crude 
a priori calculation, satellite field of view and so on. Also on the first iteration 
the initial guess of the bias is taken to be zero. The new values of 0, X and 
B obtained after the first iteration are then used as a priori 4>q , Xq and Bq for 
the second iteration. This process is repeated until the solution has converged. 
The foregoing equations can be rewritten In terms of the 3 variables S\<p^ isk 
and B. That is — 

AR = ~ + :^ AX + ~ AB A-14 

^<p-- ^X 

where 

AR = R« - R or residuals which will be minimized as the fit improves 
Aip = (p - - computed adjustment to a priori geodetic latitude 

Ax = X “ Xq = computed adjustment to longitude 

\B ■ B - Bq = computed adjustment to range rate bias (B^ taken to be zero 
on first iteration) 

In the above R Is computed on each iteration based on the current estimate 
as obtained from the previous iteration — that is. 


R 


(x - Si) X + (y “ Sj) y + (z - S3) z 




B 


A- 15 


j(x - Sj)2 + (y - + (z " S3 

as before x, y, z and x, y, z, are the satellite Earth fixed geocentric parameters 

and Sj, S^, S 3 are current beacon coordinates- The first part of the computa- 

3R 

tional problem thus becomes the evaluation of partial derivatives 


XR 


3 ’ 9X 


and which must be computed at the time of each data point collection to 
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PAGE J3 

w PocMt QUALny 


form the overdetermined set of linear equations. In this regard — 


3R 

3R 

!!? , IE 

BSj 

30 

3Sj 

30 SSj 

30 

3R 

3R 



3\ 

BSj 

3 ;^ ^ BSj 

3A. 

3R 

3B 

1,0 




aR 


p 3 


A- 16 


A-17 


A-18 


The parts of A-16».A-17, and A-18 which are not a function of the data are: 

“ V sin <p cos k A- 19 


'd<ji 


^82 

3S, 


- 1 / s in 0 sin k 


(1 - ) cos 0 


3^ 

3^ 


(v + h) cos 0 sin k - - S., 


( V + h ) cos 0 cos \ - Sj 


A-20 


A-21 


A-22 


A-23 


The parts of A- 16, A-'17, and A-18 which must be computed at each data time 
are; 


3R 

3S. 


- X + 


R(x " Sj) 


1 - 

R 


A-24 
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2L 


R(y - s,> 


-L 

R 


ii. 


_ 2 f 


R(z - Sj) 
R 


“1 


R 


A- 25 


A- 26 


In this manner equations A-16, A-17, and A-18 are evaluated at each data time 
aitd are used as the basic input to the normal equations which can be written 
as: 


‘11 t ; 


‘13 




; c, 


I ^31 


3 3 “23 


V 


A- 27 


: a 

L 


3 1 


'*3 2 ^33 i ^ 

I 




or symbolically 


where for n data points: 


A‘l) 


^ A-i C 


'12 




III 


'2 1 


A-28 


■= Ifl 
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® 13 



^31 


^2 2 



ORlGIMAi; PAGE IS 
OP POOR QUALHY 


also 


‘33 


E‘ = " 


i--= I 


^ 5k 


'23 


E 


S \. ’ “^2 

1 = 1 1 


^ , ■ 

; — AR. 

^ 5d 


4 >. 

1-1 I 


L! IX 


2 ^ 

i = l I 


^^3 - L 

i = l 

As the microcomputer steps through the data, which has first been stored 
in RAM, it pauses at each data .. le and computes the partial derivatives and 
range rate residuals and places the result in temporary storage. Each data 
point requires less than 1 second of computation time. Once the microcomputer 
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has processLni all vlata it computes the sums iiulioateil in A-liH \vhiel\ eoi'i'esptmil 


to the scalar constants Ujj throuj 5 ;h a,, auil o^, c ,, c ol A-U7, taiuulion A--7 
Is then solved for A , A\ >uul Ai^thus ohtalniuj? a solutiim for heacon locatiim 
m\d frequeuc>‘ {ranj;'e rate) bias. This takes less thai^ 10 seconds. rUe new 
values of ; , \ and B are then used as iu\ a priori for the noNt iteration attd so 
on. Althouji'h th’r i>artieular mici'ccomputer was proji'ratmued to detortnitte 
latiUule aiul longitud.e of stationary or slowlv moving beaevms itsing one or two 
IXippler data sp:u\s. the exleixsion to moving beacons and height reeiwery i.s also 
being invosiigated. If a svd’fieient niuuUiiy of iloppler data is availabU' it is in 
prineipto a simple step to n\odify the foregoing equations It' handle addtlit'ival 
m\knowns. 
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